Entanglement Entropy in Extended Quantum Systems * 
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After a brief introduction to the concept of entanglement in quantum systems, I apply these ideas 
to many-body systems and show that the von Neumann entropy is an effective way of characterising 
the entanglement between the degrees of freedom in different regions of space. Close to a quantum 
phase transition it has universal features which serve as a diagnostic of such phenomena. In the 
second part I consider the unitary time evolution of such systems following a 'quantum quench' in 
which a parameter in the hamiltonian is suddenly changed, and argue that finite regions should 
effectively thermalise at late times, after interesting transient effects. 



I. INTRODUCTION 



Entanglement is one of the most fundamental features 
of quantum mechanics, yet in some ways the most myste- 
rious. In this talk I will discuss some of the recent devel- 
opments in understanding its role in many-body systems 
extended in space, focussing on two main topics: the 
universal properties of entanglement entropy near quan- 
tum critical points, and the behaviour of entanglement, 
and, more generally, correlation functions, after what is 
termed a 'quantum quench'. 



Most of the work I will be presenting was carried out 
in collaboration with P. Calabrese P, [3] and others. 
For a more general perspective the reader is referred to 
a recent review by Amico et al. 




FIG. 1: Schematic depiction of the decomposition of the 
Hilbert space. Later these will also represent regions of d- 
dimensional space. 



II. ENTANGLEMENT ENTROPY NEAR 
QUANTUM CRITICAL POINTS 

A. Bipartite quantum entanglement 

Consider a general quantum system prepared in a pure 
state Y^), so that it has density matrix p = |vl')(vl/|. We 
suppose that the Hilbert space can be written as a direct 
product Ti = Ha €5 (see Fig. [1]). We imagine two 
observers (traditionally named Alice and Bob) such that 
Alice can make observations only in Ha (that is corre- 
sponding to linear operators of the form Oa <X) Is), and 
correspondingly for Bob. In general Alice's observations 
are entangled with those of Bob. 

One of the most useful mathematical results in under- 
standing how to quantify entanglement is that of Schmidt 
decomposition, which is based on the property of singu- 
lar value decomposition for matrices. It states that any 
pure state \^!) may be written as 



where \'ipj)A,B are orthonormal states in Ha and Hb re- 
spectively, and |cjp — 1. Moreover the Cj can be 
chosen to be real and > 0. (Note that there is only one 
sum here: for each state in Ha there is just one state in 
He-) 

One measure of the entanglement in \^!) is the entropy 
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If ci = 1 and all the rest vanish, j^*) is a product state 
and is uncntanglcd (although there may still be classical 
correlations.) If, on the other hand, all the Cj are equal, 
then S takes its maximal value, given by the logarithm 
of the smaller of the dimensions of Ha and Hb ■ For ex- 
ample, if each space is a direct product of N qubits (spin 
^ degrees of freedom) then the maximal entanglement 
entropy is A^log2. 

Equivalently, we can define the entanglement entropy 
as the von Neumann entropy Sa = —Tt^aPa^oE Pa of 
Alice's density matrix pA = Trsp. Evidently Sa = Sb = 
S. 
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Other measures of entanglement exist [j], but the en- 
tropy has several nice properties: additivity, convexity, 
basis independence. In quantum information theory 
it gives the maximum efficiency of conversion of par- 
tially entangled to maximally entangled states that Alice 
can achieve using only local operations in her part of 
the Hilbert space. From a computational point of view, 
it gives the amount of classical information required to 
specify the reduced density matrix pA- This important, 
for example, in the density matrix renormalisation group 
(DMRG). 

However, in this talk I want to use the entanglement 
entropy as a basis-independent means of characterising 
quantum phase transitions. This is particularly impor- 
tant in those cases where it is not clear what the order 
parameter is, or which correlation functions should be- 
come long-ranged at the transition. (Indeed, for topo- 
logical transitions there may be no identifiable order pa- 
rameter.) I shall consider the case when the degrees of 
freedom of the quantum system are distributed over some 
large region TZ in d-dimensional euclidean space, and the 
hamiltonian H contains only short-range interactions, for 
example a quantum spin system, or, more generally a 
UV-cutoff quantum field theory. (I shall focus solely on 
the universal properties near the phase transition: these, 
as K. Wilson and others taught us, are all encoded in the 
field theory description.) The subspace Ha will consist 
of the degrees of freedom in some large (compact) subset 
A of TZ, and we shall assume that the whole system is in 
a pure state, usually the ground state |0) of H . However, 
it will be useful also to consider the case when the whole 
system is in a thermal mixed state with p oc e~^^ . 

The main question I will address is how the entangle- 
ment entropy depends on the size and geometry of the 
region A, and on the universality class of the critical be- 
haviour. 



B. Entanglement entropy from the path integral 
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FIG. 2: The path integral in space x and imaginary time r. 
The rows and columns of the density matrix are labelled by 
the values of the fields (degrees of freedom) at t = and P 
respectively. 



z = 



FIG. 3: The partition function Z is obtained by sewing to- 
gether the top and bottom edges of the world sheet in Fig. [2] 
that is by identifying the fields and integrating them. 



analytically continue in n, and evaluate 




Sa = -Tr Pa log PA = - 



d_ 

dn 



Trpl 



For positive integer n, this is given by taking n copies of 
Fig. [4] and sewing the edges along A in a cyclic fashion, 
as shown in Fig. [5l For d = 1 (which is the case we shall 
consider here) this expresses Trp^ as a path integral on 
an n-sheeted Riemann surface, with branch points at the 
ends of the interval A (if A has several disjoint parts, 
there are several branch points.) 



"Quantum mechanics is just statistical mechanics in 
one more dimension" (M.E. Fisher) 

From an analytic perspective, the Schmidt decomposition 
is very difficult to carry out except in simple cases, and we 
employ a different route via the path integral. As Feyn- 
man taught us, the thermal density matrix cx e~^^ may 
be written as a path integral in an imaginary time inter- 
val (0,/?) (throughout we adopt units with h = 1.) This 
is illustrated in Fig. [2] The density matrix is correctly 
normalised by the partition function Z as shown in Fig. [31 
However, we need the reduced density matrix pA- This 
is found by sewing together only those parts of the upper 
and lower edges corresponding to the region B (Fig. [4]). 
In order to compute the entropy Sa = — Tr^pAlogpyi 
we use a device reminiscent of the 'replica trick' in dis- 
ordered systems: we compute Trp^ for positive integral 



1. High temperature limit 

Before proceeding further, it worth understanding how 
this formalism gives the expected answer at high temper- 
atures. Consider the case where A is an interval of length 
^ in a system of total length L, with I <^ L. Further 




FIG. 4: The reduced density matrix is obtained by sewing 
together the top and bottom edges of the world sheet in Fig. [2] 
along only the parts corresponding to region B. 



FIG. 5: A view along the edges A of n copies of Fig.|4l showing 
how they are to be sewn together to give Trp^. 



FIG. 7: An interval A of length I. When the correlation 
length ^ <^ £, the entanglement arises from regions of length 
0(5) near the boundaries between A and B. 
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FIG. 6: The world sheet for the path integral for Trp^ in the 
high temperature limit. In the central section, corresponding 
to the region A, it winds around n times. 



suppose that /3 ^ £. The path integral is then over a 
very narrow cylinder (compared to its length) as shown 
in Fig. [HI In this limit we expect the result of the path 
integral to factorisc: 

Z„«Zi(^,n/3)Zi(i 

where Zi{L,(3) is the partition function for a system of 
length L at inverse temperature f3. Thus 



exp(-n/3F^(/3)) 



where Fa{(3) is the usual Helmholtz free energy of region 
A. Differentiating with respect to n at n = 1, we then 
find 

Sa ~ KEa - Fa) . 

This shows that in this limit, the von Neumann entropy 
becomes the usual thermodynamic entropy, as expected. 



C. The critical case in d = 1 

Suppose the Id system is at a quantum critical point 
with dynamic exponent z ~ 1 (that is, linear disper- 
sion relation uj — v\k\ at low energies). Then dimen- 
sional analysis (with H = v ~ 1) implies that Fa{0) ~ 
—Trcl/GP"^. This is just the one-dimensional version of 
Stefan's law: the constant c = 1 for a single species of 
boson. We conclude that for £ ^ /3 



Tipl ^ exp 



TTC 
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The above considerations do not appear to shed much 
light on the case we wish to consider, namely zero tem- 
perature. However, in this case the limit £ ^ P is re- 
lated to the opposite one /3 ^ oo by conformal symme- 
try: the mapping z — > (/3/27r)logz converts exponential 
decay along the cylinder into power law decay at T = 0. 
Thus, at zero temperature we have 

Trp^ ^ £-{c/e)(n-l/n) ^ 



Taking the derivative we see that 

SA-{c/3)log£. 



(1) 



For a general conformal field theory (CFT), c is the cen- 
tral charge. ^ was first found in 1994 by Holzhey et al. 
0, who called it the geometric entropy. 

Note that this implies that, even at a quantum criti- 
cal point, the entropy grows only logarithmically in the 
length £, as opposed to its maximum allowed behaviour 
which is 0{£). This accounts for the success of the 
DMRG method in li = 1. This logarithmic growth is 
apparently not restricted to critical points with z — I: 
for example it also holds in random spin chains [7|. 

HJ) is just one example of a plethora of universal results 
which can be found using CFT methods. For example, 
one can explore the cross-over between finite and zero 
temperature when £ ^ (3 to find the elegant formula [l], 

Sa ^ (c/3) log {{f3/n) sinh(7r^//3)) . 



D. Finite correlation length in d = 1 

The above CFT considerations apply only at the quan- 
tum critical point. However, the behaviour of the entan- 
glement entropy close to such a point may be deduced 
from this by scaling arguments. For example, consider 
the case when A is an interval of length £ in an infinite 
system (Fig. [7]). When the correlation length ^ is finite 
and ^ £, we expect that the entanglement arises from re- 
gions of length O(^) near the boundaries between A and 
B. Since the scaling at the critical point is logarithmic 
we may therefore conjecture that, in this limit. 



5A^2x(c/6)loge. 



(2) 



and that in general there is a contribution (c/6) log^ from 
each contact point between A and B. For the case of a 
single such point, when A and B are both semi- infinite 
lines, this has been derived exactly using the corner trans- 
fer matrix [J, for a large class of integrable models. 
Recently the leading corrections to ^ have been com- 
puted These are of the form 0(e^^^/^) and they 
appear to be rather universal. 



E. Higher dimensions d > 1 

In higher dimensions, the natural conjecture for the 
entanglement entropy (given that Sa should equal Sb) 



is that it is proportional to the 'area' \dA\ of the bound- 
ary between the two regions. This has been verified in 
various models (although logarithmic factors are also pos- 
sible) and even proved rigorously 11] starting from rea- 
sonable assumptions about the behaviour of correlations. 
On dimensional grounds the coefficient of the area law 
should go as a^^'^ where is a is the short-distance cutoff, 
and therefore be non-universal. However, renormalisa- 
tion group arguments similar to those used to analyse 
the scaling behaviour of the free energy imply that there 
should be a universal term proportional to hid- 
den behind this [l|. 

It is interesting to note that these area-dependent 
terms may cancel in more complicated entanglement 
combinations. For example, if A and B are different sub- 
regions of TZ, then one may conjecture that 

Saub + Sahb — Sa — Sb 

is universal, and, at the critical point, depends only on 
the geometry and some universal constants of the critical 
theory. 



III. TIME-DEPENDENCE AFTER A 
QUANTUM QUENCH 
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B ^2t^ ^2t^ B 




2t<l 

B / B 




2t>l 

FIG. 8: Quasiparticles are emitted at t = 0, and subsequently 
move classically. Left and right-moving pairs of particles are 
entangled, and if they arrive on A and B respectively they 
cause entanglement of these two regions. 



A. The horizon effect 



In the second part of this talk I want to discuss a sub- 
ject which is intimately related to entanglement proper- 
ties of quantum many-body systems, but is somewhat 
more general. Let us suppose that we prepare a system 
at time i = in a pure state which we usually take 
to be the ground state of some translationally invariant 
hamiltonian Hq with a gap toq to the first excited state. 
For times i > we then evolve the state according to 
a different hamiltonian H (which doesn't commute with 
Hq). Note that this evolution is unitary: no dissipa- 
tion or noise. This protocol has been termed a 'quantum 
quench'. 

We may then ask how the reduced density matrix pA 
of some finite part of the system, its entropy Sa and 
correlation functions of local operators 0{x) with x £ A 
evolve. In particular, does pA reach a stationary state, 
and if so how is it characterised? 

For traditional solid state systems the assumption of 
unitary evolution usually breaks down so rapidly that 
these kind of questions have hardly been addressed in the 
past, except in a few cases of integrable spin chains p^ . 
However recent experiments on cold atoms in optical lat- 
tices have shown that it is possible to maintain coherence 
over measurable time intervals. This has prompted more 
recent work, both in integrable systems and more gener- 
ally 0, 0, Ej, [3| • These detailed considerations all lead 
to the prediction of the following rather simple physical 
effects. 



Since |^'o) has (extensively) higher energy than the 
ground state of H, it acts as a source of quasiparticles 
of H . The analysis of simple models shows that subse- 
quently they move classically. This is illustrated in Fig. [H 
for the simplest case in d = 1 when all particles travel 
at the same speed v. Entanglement between regions A 
and B (which is very small initially) arises when two en- 
tangled particles emitted from nearby points arrive on A 
and B respectively. From the upper diagram in Fig. [SI 
corresponding to early times t < £/2v, we see that the 
number of such pairs of particles, and therefore the de- 
gree of entanglement, increases linearly with t. One the 
other hand for t > £/2v is saturates at a value propor- 
tional to i. This is what is found in explicit calculations. 
The coefficient of the linear term depends on the initial 
state. 



1. Correlation functions 

In general it is found that one-point functions {0{x, t)) 
of local operators (for example spins in a quantum spin 
chain where |^o) breaks the spin-reversal symmetry but 
H and its ground state do not), decay exponentially fast 
towards their values in the ground state of H. The 
physical reason for this is illustrated in Fig. [HI On 
the other hand the behaviour of a two-point function 
{0{xi,t)0{x2,t)) is governed by the horizon effect. Up to 
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FIG. 9: A local observable, for example a quantum spin, is 
subject to incoherent radiation of quasiparticles emitted from 
different spatial points. Its value therefore flips according to 
a Poisson process, and its expectation value relaxes exponen- 
tially. 




FIG. 10: Two examples of dispersion relations. In each case, 
the quasiparticles move with their group velocity given by the 
slope of this curve. The fastest moving particles are respon- 
sible for the horizon effect, and the slowest particles for the 
approach to the limiting behaviour at late times. 



time t ~ \x2 —Xi\/2v its connected part does not change 
from its initial form which is very short-ranged. For wider 
separations {0{xi,t)0{x2,t)) ^ {0{t))^, which decays 
exponentially in time. At the time i ~ |x2 — xi\/2v the 
two points fall inside the horizon and the connected cor- 
relation function becomes non-zero. In the case when 
all quasiparticles have the same speed, the full correla- 
tion function then becomes time- independent, and, since 
it was previously decaying exponentially with t, now ex- 
hibits exponential decay in the separation \x2 — xi\. 



2. General dispersion relation 

These physical considerations allow us to understand 
the behaviour in more realistic situations where lattice 
effects or a gap modify the dispersion relation. Exam- 
ples are shown in Fig. [TOl The fact that the quasiparti- 
cles, once emitted, travel classically means that it is their 
group velocity which is important. Since this can van- 
ish, for example at the zone boundary, this is responsible 
for a very slow (power law + oscillations) approach to 
the asymptotic limit at large t, which is actually seen in 
exact calculations and numerical work. 



B. Thermalisation at late times 

We have argued that calculations in various solvable 
models suggest that correlation functions for x £ A, and 
by implication the reduced density matrix pA, become 
stationary at late times. But what is their form? In these 
cases it turns out that, at least in the case where the gap 
mo in the spectrum of Hq is large, pA has a thermal 
form oc e^^""^ , where (3~^ is an effective temperature 
dependent on mo among other things. 

In order to understand this, first consider a very sim- 
ple example: a simple harmonic oscillator quenched from 
frequency ujq (and initially in the ground state), to fre- 
quency uj. It is a simple calculation to evaluate the over- 
lap between the ground state |5'o) and a typical eigen- 
state of H with energy E. For o^o ^ w we find 

(*o|S) cxexp(-/3effS/2),, 

where (ScS ^ 4/wo, so the matrix elements of the time- 
dependent density matrix have the form 

{E\pit)\E') ^ e-/3c«(B+£')/2g.(B~B')* . 

Note that since the energy differences are all multiples 
of Lu, p{t) does not tend towards a stationary value - it 
oscillates as expected! 

However in an extensive system, which we can consider 
as a set of oscillators, one for each quasiparticle mode fc, 
if we consider a finite region A of size ^ t/w as i — *■ oo, we 
need to integrate over all the fc-modes. If the dispersion 
relation has the form LOk = m + 0(fc^) the different modes 
destructively interfere except at fc = 0. The result is 
that correlation functions, and pA{t)^ become stationary 
as i — > oo as if they were at finite temperature. In fact 
a more careful calculation, valid not just in the large mo 
limit, gives 

PcS — (4/m) tanh^"'"(m/mo) . 

Similarly, we find that the extensive part of Sa saturates 
at a value equal to the thermodynamic entropy at this 
effective temperature. 

Although the above simple argument relied on a non- 
zero gap m > 0, in fact (as long as interactions are 
present) the above result appears to hold (for example 
in a CFT in 1+1 dimensions) even in the gapless case. 

It should be stressed that this effective thermalisation 
occurs despite the fact that the system as a whole remains 
in a pure state, and there is no ergodicity or coupling to 
a heat bath - the effect arises solely as a consequence of 
quantum interference and entanglement. 

Summary 

In this talk I have argued (a) that entanglement en- 
tropy provides a useful order-parameter independent di- 
agnostic of quantum phase transitions, with many uni- 
versal features, and (b) that after a quantum quench. 
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there are interesting transient phenomena hke the hori- 
zon effect, and that at late times finite regions should 
behave as though they are in thermal equilibrium. 

There are many open questions, particularly in regard 
to part (b), and whether it holds only for theories which 
admit a quasiparticle picture. In addition, these kinds of 
calculations need to be repeated for more realistic mod- 



els, and also for quenches through the critical point into 
an ordered phase, where the relation to later-time coars- 
ening effects needs to be understood. 
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